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Abstract. We give new constraints on small-scale non-Gaussianity of primordial curvature 
perturbations by the use of anisotropies in acoustic reheating. Mixing of local thermal or local 
kinetic equilibrium systems with different temperatures yields a locally averaged temperature 
rise, which is proportional to the sqnare of temperature perturbations damping in the photon 
diffusion scale. Such secondary temperature perturbations are indistinguishable from the 
standard temperature perturbations linearly coming from primordial curvature perturbations 
and hence should be subdominant compared to the standard ones. We show that small- 
scale higher order correlation functions (connected non-Gaussian and disconnected Gaussian 
parts) of primordial curvature perturbations can be probed by investigating auto power 
spectrum of the generated secondary perturbations and the cross power spectrum with the 
standard perturbations. This is simply because these power spectra come from higher order 
correlation functions of primordial curvature perturbations with non-linear parameters such 
as /nl and tnl since secondary temperature perturbations are second order effects. Thus, 
the observational results l{l + ~ 6 x 10“^*^ at large scales give a robust and universal 

upper bound on small-scale non-Gaussianities of primordial curvature perturbations. 

Keywords: scale-dependent non-Gaussianity, primordial bispectrum, primordial trispec¬ 
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1 Introduction 

The 21st century is the age of precise observational cosmology. It is possible to observation- 
ally test theoretical models for the early universe by analysing anisotropies of the Cosmic 
Microwave Background (CMB) [1, 2]. In fact, we already know that primordial curvature 
perturbations do exist, and their spectrum is almost scale-invariant. While the statistical 
property of primordial perturbations is well approximated by Gaussian statistics, the devia¬ 
tion from the exact Gaussian statistics [3] will deliver us rich information on the primordial 
universe [4], which is one of the next targets in cosmology. On the other hand, it should 
be noticed that the scales probed by the CMB anisotropies are just 6 e-foldings of the last 
60 e-folds during inflation and that we know very little about smaller-scale perturbations 
on k Mpc > 0.1 because most of the primordial density fluctuations dissipate due to Silk 
damping. One powerful tool to see the fluctuations on small scales is CMB spectral distor¬ 
tions [5-13]. Deviations from the ideal blackbody spectrum are induced in the early universe 
after around z ~ 2 x 10®, and such deviations are characterized in terms of several parame¬ 
ters such as the Compton y parameter and the chemical potential y. The next observational 
projects, such as PIXIE and PRISM [14, 15], following COBE/FIRAS are proposed [16-18], 
and it may be possible to see the nature of primordial fluctuations up to another 11 e-foldings 
in the near future. Hence, 17 e-foldings will be within the scope of our current technology, 
but more than one half of total e-foldings necessary to solve the initial condition problems 
remain obscure. 

Recently, acoustic reheating was proposed to investigate primordial curvature perturba¬ 
tions on extremely small scales [19, 20]. Below the diffusion scale, the inhomogeneities are 
erased due to Silk damping so that the universe is heated by the conversion of the dissipated 
energy of photon [21]. Such temperature rises are quadratic in temperature fluctuations over 
the diffusion scale. In the previous works, the total average (homogeneous part) of the tem¬ 
perature rise was discussed, and they derived the constraints on small-scale power spectrum 
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of primordial curvature perturbations by comparing the helium mass fraction Yp (or pho¬ 
ton baryon ratio tjb) during Big Bang nucleosynthesis with that at the last scattering of the 
CMB. Such investigations are new and interesting, but the scales probed by these approaches 
are limited up to /cMpc ~ 0(10®) since the heating in much smaller scales (at much earlier 
epochs) changes both and so that we know very little about the scale except 

the constraints that Vti ^ 0.3 for k ~ 10^^“^®Mpc“^ [19]. However, the temperature rise in¬ 
duced by acoustic reheating is only homogenized over the diffusion scale. Then, it fluctuates 
spatially in general and potentially includes all of the histories of erased temperature pertur¬ 
bations from the beginning of the hot universe, that is, if large inhomogeneities would exist at 
small scales, the observed temperature anisotropies can be changed drastically. In this paper, 
we shall focus on such inhomogeneities of acoustic reheating. Although one cannot pick up 
and isolate the temperature perturbations due to acoustic reheating from the observed tem¬ 
perature anisotropies, given the fact that the linear theory prediction well matches with the 
CMB observations, the non-linear corrections should not be beyond the observational Cf"^, 
i.e. the order of 10“^*^ at large scales. This puts robust and universal constraints on higher 
order correlation functions including non-Gaussianities of primordial curvature perturbations 
since the non-linear temperature corrections are of the second order of primordial curvature 
perturbations. Our approach is closely related with a method to investigate anisotropies in 
CMB distortions which was originally proposed by Pajer and Zaldarriaga since the physics 
of acoustic reheating is essentially the same [22]. 

We organize this paper as follows. In section 2, we summarize the details of acoustic 
reheating, particularly paying attention to its perturbations. Section 3 and 4 are devoted to 
the calculations of angular power spectrum by the use of Sachs-Wolfe approximation, and 
we obtain the constraints on the non-Gaussianities of primordial curvature perturbations in 
section 5. We give conclusions and discussions in the final section, and scale-dependence of 
non-Gaussianities is briefly mentioned in the appendix. 

2 Inhomogeneities in acoustic reheating 

The local photon temperature in the early universe fluctuates due to the primordial curvature 
perturbations. As long as the Compton and the double-Compton processes are efficient, 
which holds for z > 2 x 10® [23, 24], it is locally in thermal equilibrium, and its spectrum 
obeys the Planck distribution characterized only by the local temperature T{r],x,h), where 
r], X and fi are the conformal time, the space coordinate and the unit vector of the photon 
momentum respectively. This local temperature can be divided into the homogeneous part 
and the inhomogeneous one. 


r(77,®,h) = r(7/)(l -k 0(r?,a;,h)). (2.1) 

Here, the acoustic reheating is not yet taken into account, and hence T decreases in propor¬ 
tional to a{rj)~^ (o(^) : the scale factor) in the expanding Universe as long as the number 
of the effective degree of freedom of radiations is unchanged, and the temperature perturba¬ 
tions are proportional to the primordial curvature perturbations at linear order. Once the 
acoustic reheating happens, small-scale inhomogeneities are erased due to Silk damping, and 
the effective temperature is raised because of the energy conservation. Now, let us evaluate 
the temperature rise for a given perturbed system of photons, assuming acoustic reheating 
happens instantaneously. Keeping in mind the relation between the temperature and the 
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energy density of photons p, T = {^pjaB)^^^ with the numerical constant as = jlh, the 
dimensionless temperature rise A is given by 


A(a;) 



( 2 . 2 ) 


where {X)x represents a spacial average defined as 

{X)x = [ d^x' Wrj,{x')X{x + x') 


(2.3) 


with Wrj, being a window function with a radius vt around x. This spacial averaging 

procedure effectively represents the (instantaneous) acoustic reheating, in which we assume 

that damping and thermalization happen instantaneously within the same box. The term 
1 

{p)x represents the average temperature when the energy density of photons around x would 
be damped perfectly and thermalized homogeneously within the radius r^. Another term 
{p^)x represents the average temperature when the (fluctuating) photon temperature would 
be homogenized within the radius rT- From the energy conservation, this difference yields the 
temperature rise. One can easily see that this effect comes from nonlinear relation between 
the energy density and the temperature of photons, and is second order in perturbations of 
the photon temperature. In fact, inserting the definition of the photon perturbation (2.1) 
into (2.2) yields A(a;) = 3(0?)a;/2. Here 0j represents the temperature perturbation before 
acoustic reheating happens. Realistic acoustic reheating proceeds through two processes, 
damping of the photon perturbations due to Silk damping followed by (homogeneously) 
thermalization through the Compton and the double-Compton processes. Therefore, the 
actual temperature rise for a conformal time rj can be written as 


= -(0| - Q {v))x- 


(2.4) 


Here, we separate two processes explicitly. 0? — 0^(7?) represents only the damping effect and 
includes the time evolution of damping scales. The spacial average in this equation represents 
only (homogeneous) thermalization effect and does not include the damping effect. Then, 
this equation can be recast into the evolution equation of A as 


dA 

dr] 


3 d{e^) 

2 dp 


(2.5) 


Strictly speaking, thermalization (diffusion) radius also depends on the conformal time. 
However, as long as we are interested in the final temperature rise, we have only to take 
to be the largest thermalization radius during acoustic reheating. Therefore, we fix rx to be 
such a radius here and hereafter. 

During 5 x 10^ < z < 2 x 10®, the double Compton process is negligible while the 
Compton process is still efficient. Then, the system is locally in kinetic equilibrium, but its 
spectrum obeys the Bose-Einstein distribution with a nonzero chemical potential p [23-25], 
which modifies the expression of the temperature rise A for a given perturbed system of 
photons, assuming acoustic reheating happens perfectly and instantaneously, as 


A = 


{p)x - {p*)x _ 45C(3) 
f 27r4 




( 2 . 6 ) 
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with C being a Riemann zeta function. Subtraction of the chemical potential means that the 
temperature of the Bose-Einstein system is not the forth root of the energy density. Then, the 
relation between the temperature rise, and the photon perturbation is also modified as [13] 

== - 2 . 277 '^ 

dr] dr] 

Thus, in order to estimate the amount of the acoustic reheating, we have only to evaluate how 
much the temperature perturbations originated from the primordial curvature perturbations 
decay due to Silk damping. Once the acoustic reheating is taken into account, the local 
temperature is modified from (2.1) into 

r(ry, X, h) = T{r]) [l + 0(ry, n) + A(77, x)]. (2.8) 

It should be noticed that A consists of only a monopole component because of the ther- 
malization effects and that it will include the homogeneous part for the entire Universe in 
general, which can reheat the whole Universe. 


(2.7) 


The time evolution of the linear temperature perturbations in the conformal Newtonian 
gauge is described by the Boltzmann equation [26-29], 


d^ 

dr] 


0 


(f) — {k ■ n)'ij) + f 


0-00+2^2(fc 


n)(02 + Q2 + ^0) ~ 


(2.9) 


where r is the optical depth, and fc • n is a cosine between unit vectors of the photon 
momentum and the Fourier momentum. 0ver-dots represent partial derivatives with respect 
to the conformal time. Each quantity is Fourier transformed, and 0 is also expanded by 
Legendre polynomials Pi as Q = X]z(“*)^(2^ + ' h.)0;. 0^ is the polarization, and its 

Legendre coefficients are given in the same manner. The definitions of the velocity potential v, 
gravitational potential ^|J and the curvature perturbation (p are based on [28] . The third term 
on the r.h.s. is negligible on large scales. Then, the monopole and the dipole components of 
temperature perturbations are enhanced due to the Sachs-Wolfe effects just after the horizon 
entry. 0nce (p and ^p decay, the monopole and the dipole oscillate in the inverted phase. 
The leading term proportional to r comes from the quadrupole component, and emerging of 
such nonzero anisotropic stress induces the diffusion damping. Below the diffusion scale, the 
Boltzmann hierarchies are solved approximately [30], 

1 P k‘^\ 

01 ~ sm{krs) exp (2.10) 

where TZk is primordial curvature perturbations on comoving slice, kp^^ corresponds to the 
diffusion scale, and is the sound horizon. In the quasi tight coupling regime, the following 
relations approximately hold, 0^ + Qq = §02) ©2 = 8A:0i/(15'f) and drjkp^ = —8/(45f). 
As pointed in [19], the diffusion scale should be divided into at least 5 regions depending 
on redshift dependence of /cd and the chemical potential. Taking the beginning of the hot 
universe to be zq = 10^*^, in which anisotropic shear comes from relativistic species such as 
7, and Z bosons until zi ~ 1.5 x lO^'^, damping scale is written as ko — 5.5 x 10®(1 + 
_j,)0.5iMpc“^. 0nce the weak bosons acquire their masses, neutrino is a dominant component 
to smooth inhomogeneities before neutrino decoupling around Z 2 ~ 5.5 x 10®, and the diffusion 
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scale is given by ko — 5.Ox 10 “^^(l+ 2 ;)^'^Mpc“^. Since neutrino can smooth inhomogeneities 
over the whole horizon scale just before the decoupling, the diffusion scale is almost the same 
with that at decoupling. Then, the diffusion scale is almost constant until the photon diffusion 
scale goes beyond ko ~ lO^Mpc”^ around ~ 8.4 x 10®. zs ~ 5 x 10^ is the lower bound 
since kinetic equilibrium is not established anymore. Thus, during < z < z^, the damping 
scale is induced by photon diffusion and is written as ko — 4:.l x 10“®(1 + 2 ;)®/^Mpc“^. 
Though the diffusion scale is unchanged, we need to distinguish a chemical potential era 
from a blackbody era, whose transition happens around 2:4 ~ 2 x 10®. Combining (2.5), 
(2.9), and (2.10), the temperature rise in the Newton gauge due to the acoustic reheating at 
the last scattering conformal time r/* is evaluated as 


A(r/*,a;) 



d^k2 
(27r)3 (27r)3 


e*('“i+^^)'®VV’rr(fc)5T2(fci • n)P2{k2 ■ fi){sm{kirs)sm{k2rs))p 


5 

X ^ ^ 0!n 
n=l 



kf + ^2 

khiz) 


^n — 1 


( 2 . 11 ) 


where ■ ■)p is periodic average over the duration longer than the oscillation period. Zn 
and an is given by {zq, • • • , 25) = (10®®, 1.5 x 10^^, 5.5 x 10®, 8.4 x 10®, 2 x 10®, 5 x 10^) and 
(ai,--- , 04 , 05 ) = (3/2,-- - ,3/2,2.28), respectively. In (2.11), we drop the dipole contribu¬ 
tion proportional to 0 i 0 iPi(fci)Pi(fc 2 ) because its periodic average vanishes, and omitted 
the higher order multipoles proportional to PiPi' with I ^ I', which vanish thanks to the 
orthogonality of the Legendre polynomials when we estimate the angular power spectrum 
later. Coarse graining is operated through the window function Wr^i^k), which is the Fourier 
transformation of Wrj,{x) defined in (2.3). The typical scale of such homogenization is given 
by the maximum scale of thermalization around z ~ 5 x 10^. Then, inhomogeneities in 
acoustic reheating are erased on scales below O(O.lMpc). 

In Fourier space, (2.11) can be written as 


A(r/*,fc) =2^ On J j ^^^5P2{ki-n)P2{k2-n){2Trf6^^\k-ki-k2) 


Wrj, (k) {sm{kirs) sm{k2rs))^ 


exp 


khiz) 


n Zn-l 


TlkiT^k2 


( 2 . 12 ) 


Here, it should be noticed that the temperature rise due to acoustic reheating is smoothed 
over the diffusion radius so that the contribution for k ^ is significantly suppressed. 
However, this fact implies that only the sum of ki and k 2 must be small and that ki and 
k 2 themselves can be large because the acoustic reheating is the second order effect. This is 
the essential reason why the perturbations of acoustic reheating can probe the information 
of primordial curvature perturbations on small scales. 

Before the end of this section, we shall make a comment on the other 2nd order effects. 
In the above 0 is assumed to be given by the solution of the linearized Boltzmann equation, 
(2.9). As studied in [31-35], the non-linearity of the Boltzmann and the Einstein equations 
also induces non-linear temperature fluctuations of the CMB, namely the non-linear part of 
0. While such 2nd order effects are so far evaluated around and after the recombination, 
we investigate a part of such 2 nd order effects generated well before the recombination in 
this paper. Then the effect estimated in this paper will be separately treated from the 
standard 2nd order effects simply because the physical origin is different. Of course, in order 
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to compare with the observed temperature anisotropies of the CMB, one should take into 
account the both effects, which are indistinguishable in fact. Here, to emphasize and note 
the importance of 2nd order effects well before the recombination, we shall especially focus 
on this part of specific 2nd order effects. 

3 Angular power spectrum 

Temperature anisotropies which we observe are decomposed using the spherical Harmonics 
Yim, and their coefficients are defined as 

= j dfiYi*^{h)Q{x = 0, n), (3.1) 

aL = J dnYi*^{n)A{x = 0,n), (3.2) 

where we have set the coordinate of an observer to the origin without loss of generality. 0 
and A above are evolved from the last scattering surface and hence A has h dependence 
as well in contrast to (2.11). The coefficient of the total temperature perturbations is given 
^Tm ~ ^fm + Then, the CMB temperature angular correlation is written as = 
C®® + 2(7®^ + where we have defined 

(3-3) 

m=—l 

and (• • •) represents the ensemble average. From the observations of the CMB anisotropies, 
at large scales are known to be the order of 0(10“^*^). Therefore, unless each contribu¬ 
tion cancels out accidentally, it should be smaller than ©(lO”^*^), that is, C®^ < 0(10“^*^) 
and < 0(10“^^). These conditions yield new constraints on the primordial curvature 
perturbations at small scales. This is the central topic of this paper. 

Using the Fourier modes, we can express (3.1) and (3.2) as 

/ d^k 

j^Y,Uk)Ti{k,7jo)e{7i*,k), (3.4) 

/ d^k 

j^yi*mik)Ti{k,7jo)A{rj^,k), (3.5) 

where r]Q is the present conformal time. The transfer function is given by 7; — [0o(??*) + 
{'n*)]ji[k{do ~ V*)] — ~^ji{kr]Q)/5. We can use the same transfer function for 0 and A be¬ 
cause they are indistinguishable. Using a relation between spherical harmonics and Legendre 
polynomials 


^ 9/11 

YC^ik,)YUk2) = • ^ 2 ), 


(3.6) 


(3.3) is reduced to 

= ^ J jf{kr]o)PxY{k) = ^ J dink jf{kr]o)VxY{k). (3.7) 
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Here Pxy and Vxy are dimensionful and dimensionless power spectra for 0^ = 0(r/*,fc) 
and Af^ = A(? 7 *, k) and are defined as 

{XkY^,} = (27r)35(3)(fc _ k')PxY{k), VxY{k) = ^PxY{k), (3.8) 

where X, Y is either 0^ or A/,. 


4 Power spectrum 


In this section, we estimate the cross- and auto-correlations, namely (0fcAfc/) and (A^Afe/), 
to calculate the non-linear corrections to the temperature angular power spectrum such as 
and in the next section. 


4.1 0-A cross power spectrum 

First of all, (2.12) yields the following cross correlation function: 




E' 

n=l 


d^ki j" d?k2 

(27r)3 J (27r)3 


{27T)^6^^\k' + ki + k2)bP2{ki ■ h)P2{k2 • n) 


X Wrj,{k){sm{kirs)sm{k2rs))f 


exp 


kf + ^2 

k%{z) 


^n — 1 


{'JZk'Pki'Pk2 


(4.1) 


In this paper, we concentrate on local type non-Gaussianities because large- and small- 
scale perturbations must be correlated to probe small scale inhomogeneities through acoustic 
reheating on small scales, by using the CMB anisotropies on large scales. Then, we shall 
replace the bispectrum of primordial curvature perturbations by the scale dependent 
non-linear parameter /nl with the aid of (A.l) and (A.2) in the Appendix: 


4 ^ 

P0a(/c) ^->VrT(fc)x”" ^Vniko)^ar, 

n=l 


dlnki fNL{ki,ki,k)x'^‘’ ^ exp ( - 


2kl 

khip 


Zn-1 


(4.2) 


where we have assumed <C fei because we are interested in the corrections coming from 
small-scale perturbations to the large-scale CMB anisotropies. Therefore, this expression is 
exact only in the limit of k —)• 0. For large k modes, the additional suppression factor 
exp {—k‘^/kjj) appears. Here and hereafter, we assumed n* < 4 otherwise stated, x and xi 
are defined hy x = k/ko and xi = ki/ko respectively where ko represents a wave number at 
the pivot scale. Below we shall consider two simple but interesting types of scale-dependent 
non-Gaussianities, which are encoded in the scale dependence of /nl parameters as discussed 
in the appendix: the power law type and the top hat type. 

Let us consider the power law type given in (A.4) and (A.5). For this type, the scale 
dependence is defined by geometric and arithmetic averaged wavenumber, whose cross power 
spectra are given by 




^fSL^(k)WrAk)x'“ ‘PK(t(])r (e°'''‘) E 


O^r. 


n=l 


kpjz) 

2kl 


.g/a- 


(4.3) 
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where GjA indicates geometric average and arithmetic average respectively, and r(x) rep¬ 
resents the gamma function. For the former case, /nl(^) = ~ + 

2nfl3 — l)/2 while, for the latter case, /nl(^) = /n^^( 2/3)"’-^ and = {ug + nj — l)/2. 
Bluer spectrum index induces larger temperature rise, and this power spectrum is almost 
scale-invariant on large scales because their spectral indices are given byns-|-n//3 — 1 and 
Us — I respectively. 

Next, let us consider the top hat type. Substituting the top hat type /nl defined by 
(A. 10) into (4.2), we obtain 


^tophat(^) 


= ^yVrj.{k)'P^{ko)x'^‘’ 

n=l 


^CMB 

Jnl 



+/nl 


.ns -1 

i 


El 

2 



— X 


ris — l 

f 


El 

2 


(3-na) 



^n — 1 


Zn 


(4.4) 


where E„(x) is the exponential integral function defined as 


E 


n 



dt 



(4.5) 


The first term represents the contribution coming from the scale-invariant part Ey 

taking into the present constraint on the (scale-independent) local type /nl, this term is 
easily shown to be negligible. This power spectrum is also almost scale-invariant on large 
scales because its spectral index is given by — 1, and we can see that the amplitudes are 
proportional to the logarithmic width of the top hat given by logki/kf since Ei{x) ~ logx. 
This implies that even if we take k~^ as small as the horizon at the beginning of the universe, 
the total contributions to acoustic reheating are not significant but at most the logarithm of 
the fraction kj and ki. 


4.2 A-A power spectrum 

In this subsection, the power spectrum of A is computed, and there are 3 types of con¬ 
tributions: two non-Gaussian connected parts and Gaussian disconnected part. Therefore, 
VAA{k) is nonzero even if the primordial curvature perturbations are purely Gaussian. Using 
(2.12), the auto-correlation function (AfcA)(,,) is written as 



X 

X 

X 

X 


5P2{ki • fi)P2{k2 ■ n)5P2{ks ■ n)P2{ki ■ fi) 
(2^f6^^\k -ki- k2){2TTf5^^\k' + fcs + k^) 

WrT{k){sm{kirs)sm{k 2 rs))p exp — hrT^ 

r / _j_ jij.2 

Wrj.{k'){sm{k 3 rs)sm{k^rs))p exp —rUTV 

L V kfy{z) 


Zn — 1 

Zn 

Zm — 1 

Zvn 


X {'Pki'Pk2'^k3'Pk4) 


(4.6) 
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Ensemble average of the fourfold product of TZ can be divided into connected and disconnected 
parts, and the former can be parametrized as 


{TZkJZk,2TZkzTZk^) c ={2tt)^ + ^2 + % + ^ 4 ) 


X 


'rNL(A:i, k 2 , h, ki){Pn{ki2)Pn{ki)Pn{k3) + llperms.) 


+ ^ 2 , fcs, k4){PTi{ki)Pn{k2)Pn{kz) + Sperms.) 


(4.7) 


On the other hand, using Wick’s theorem, the latter is reduced to 


{nk,nk2nk,nk,)dc =(27r)35(3)(fc^ + k2){27:f6^^\k3 + ki)PTi{ki)PTi{k3) 

+ {2nf6^^\ki - k3){27rf6^^\k2 - fc4)P7^(^l)^’7^(^2) 

+ {27Tf6^^Hki - fc4)(27r)3<5(3)(fc2 - k3)Pn{ki)Pn{k2) 


={2T:f5^^\ki + k2-k3-ki) 


{27:f5(^\k3 + k^)Pn{ki)Pn{k3) 


+ {2T,f5^^\k2 - k^)Pn{ki)Pn{k2) + {2TTf5^^\k2 - k3)Pn{ki)Pn{k2) 

(4.8) 


Disconnected part 

Collecting disconnected contributions, (4.6) and (4.8) yield 


x^l dlnki 



■2^71 — 1 

■2^71 


2 


(4.9) 


where we have omitted the term proportional to 6^^\k) because such a term does not con¬ 
tribute to except the monopole, which should be included in offsets. It should be again 
noticed that we have assumed <C fei in the same way as the 0-A cross power spectrum. 
Therefore, this expression is exact only in the limit of A; —0. For large k modes, the addi¬ 
tional suppression factor exp {—k'^/kj^) appears. After the momentum integration it yields 


^AA(fc) W' 


rjp 


ik)Pl{ko) 2^ 


- 1 


X r 


2n, — 5 


a. 


K-O 


kn(z5) 


(4.10) 


where we have used k^izi) > koizj) for Zi > Zj. If we assume the scale-invariant power 
spectrum with Vn{ko = 0.05Mpc“^) ~ 2.2 x 10“®, (4.10) yields 


V^^{k ^ 0) 61.7 X 


koiz^) 


3 

Pn(^o)- 


(4.11) 


This implies that the spectral index of large scales is 4, and the pivot scale is k^^z^) S> 

ko- Therefore, very little contributions to the COBE scale temperature perturbations are 
expected. 
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Connected part 

Among connected contributions, the most important contribution comes from the term with 
tnl; which is written as 


^aa(^) / dln/ci / d\iik2 mh{ki,ki,k2,k2)V'ji{ki)Vn{k2) 


E 


OLr, 


exp 


2kl 

kl{z] 


^n — 1 


E 


Or 


exp 


2k‘l 

k^iz) 


^m—1 


(4.12) 


where we have used (4.6) and (4.7), and have assumed A: <C fei, ^2 in the same way. We shall 
particularly focus on the geometric averaged type for simplicity because the expression for 
the arithmetic averaged one is a bit redundancy. For the geometric averaged type, using 
(A.9), the power spectrum reduces to 



(4.13) 


where e = (n^ + nr/2 — l)/2. The scale dependence is similar to that of and is almost 
scale invariant on large scales. On the other hand, in the case with the top hat type non- 
Gaussianity given in (A. 10), (4.12) yields 


v2T{k) =vl{ko)wl{k)x 


CMB-p2 I IT-s — ^ 


E' 

k n=l 


ns — l-\ ^n—1 ' 

2k Q 


+ 2rNL F 


n.c - 1 


E' 

. n=l 


X E 

\m= 

+^nl 


(Xr 


m=l 

5 

(Xr, 

m=l 


ns —1 t7 ^ 

Xz 11/3 —ns 


ris-l-p 
Xz ihS—ns 


kljz) 

2kl 

J±_ 

2k%{z[ 


2kl,{z) 


^n—1 ' 


ns —1 t7 ^ 

X j; H/ 3 — ns 
J 2 


X f Hj 3—ns 
J 2 


k} 


2kl,{z) 


1-2 


^n — 1 '' 


-2^n —1'' 


2k‘j-){z 


(4.14) 


In this case, the scale dependence of is also similar to that of and is almost 

scale invariant on large scales. Finally we shall also show the another contribution from the 
term with q'nl for the sake of completeness: 


-^r4(^o)rV’r^(fe) / dln^i / dlnfc2 ^ 


E 


(Xr. 


exp 


2ki 

khiz. 


^n — 1 


E 


(Xr, 


exp 


2/cl 
kl){z) 


^m — l 


(4.15) 


Note that the spectral index of this power spectrum is 4 as in the case with the disconnected 
part and extremely blue on large scales. 
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5 Non-linear corrections to the temperature anisotropies and constraints 
on non-Gaussianities 


Unfortunately, the scale-independent non-Gaussianities of primordial curvature perturbations 
which are now strongly constrained by the observations of the CMB anisotropies have little 
contribution to non-linear temperature corrections. However, the non-Gaussianities might 
not be (even approximately) scale-invariant over a wide range of scales, which implies their 
magnitude on small scales can be large (or small). As shown thus far, the non-linear correc¬ 
tions to the temperature anisotropies from acoustic reheating are sensitive to the local type 
scale-dependent non-Gaussianities. Then, in this section, we concretely estimate the GMB 
anisotropies coming from such temperature perturbations and give constraints on small-scale 
non-Gaussianity of primordial curvature perturbations. Our purpose is not full analysis but 
to obtain the Sachs-Wolfe plateau for the low I’s. Therefore, we can safely set yVrj.{k) = 1 
because krx 'C 1 and omit the exponential suppression factor due to the photon diffusion 
hereafter (except the disconnected and the (7 nl cases as explained later) since the spherical 
Bessel functions only project small k modes onto low I plateau. 

The angular cross power spectra between 0 and A in the case of power law type non- 
Gaussianity are obtained from (3.7) and (4.3) as 




=^Y^/Nr^l(A:o)r(6«) a. 


n=l 


U2 \ 

2kl 


Gn Zn-l 


G{l,ns -t- n//3). 


C 


0A,A 97r2 /2\"^ 


125 V3 


/CMBp|(A:o)r(e^)^a, 


n=l 


2k Q 


A-\ Zfi—\ 


0(1,Us), (5.1) 


where we define 


on r(/ + t-^)r( 3 -n) 


(5.2) 


Here we have used the expressions of the power spectra only for small k modes, which is 
justified for our purpose because the spherical Bessel functions only project small k modes 
onto low I plateau. For the top hat case, (3.7) and (4.4) yield 


c; 


©A,tophat 

I 


-Y^h-LV^{ko)g{l,ns) 


X ^ ^ C^r, 
n=l 




2kf 




2k'j 


Zn-1 


(5.3) 


Angular auto power spectra of A with different scale-dependent non-Gaussianities are also 
calculated from (3.7), (4.13) and (4.14). The formulas are summarized as 




Ci''^^^T§}^^VUko)Q{l,ns)[T{e)Y,o^r^ 


n=l 


u2 \ e- 
Hd 


2kl 


Zn-l ' 


(5.4) 
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and 


^AA,tophat 
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Zn-1 
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.ns-1 
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3 —ns 
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— X 


ns —1 
/ 


E 3—n, 
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+'rNL 




El 

2 



— X 


ris-l 

/ 


El 


(3-ns) 



-2n— 1 


Zn 


2 ' 


(5.5) 


where we have omitted the first term in (4.14) because it is negligible. 

Now that all the formulae are derived, let us investigate the physical implications of the 
results. As we have already stated before, given the fact that the linear theory prediction 
well explains the observations of the CMB anisotropies, the corrections obtained above should 
not exceed the observed values and hence must be subdominant. From this condition, we 
obtain the constraints on small-scale non-Gaussianities. In Fig. 1, we show the constraints 
on the non-Gaussianities with the power law type of the geometrical average at pivot scale 
ko = 0.05Mpc“^ as a function of each spectrum index with fixing n* = 1. Horizontal axes 
represent the spectrum indexes of the non-Gaussianities, and the vertical axes represent the 
logarithms of the non-linear parameters. The contours of the /Cff^/Cf®) are 

inserted at the intervals of 10^ in the left (right) panel, and red dashed lines correspond to 
the lines along which the non-linear corrections are comparable with Cf®. Then, only the 
region below the red line is allowed. 

In fig. 2, we show the constraints on the top hat type non-Gaussianities in the case 
with Us = 1. For simplicity, we also set tnl = 0. Instead of giving the constrairhs on /nl 
and Tnl, we give the constraints directly on the magnitudes of bispectrum B-ji = /nlT’.^(A:o) 
and trispectrum Tti = T^\J^^{ko). This is because /nl and tnl are normalized not by the 
(unknown) power spectrum at the corresponding scale but by that at the large (CMB) scale, 
that is, P'ji{kQ). The figure shows the constraints on bispectrum B-fi and trispectrum Tti. 
The contours of (2(7^^ -|- are inserted at the interval of 10 in the case with 

logi(^ ki/kf = 4, and a red dashed line corresponds to the line along which the non-linear 
correction is equal to Then, only the region below the red line is allowed. 

The disconnected part and the non-linear term proportional to ^'nl do not contribute 
significantly. This is because such configurations do not pick up a pair of perturbations which 
are confined in the diffusion scale. The non-linear corrections to angular power spectra coming 
from these parts become 


^AA,dc/g 


367r r 

J 

SOvr 

■25(27r2 


d^k 

(27r)’ 


jf{kr]o)PA^^ik) 


P^/^(fc^0) / dkejfikrio)WrAk)- 


(5.6) 


Note that, in these cases, the use of the power spectra only for small k modes is not justified 
because their spectra indexes are 4 and extremely blue, as already pointed out in (4.10) and 
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Figure 1. The constraints on plane (left) and plane (right). The interval of the 

contours of the 2C^"^/C'^®(left) and (right) are 10®, and red dashed lines correspond to 

the lines of 2(7^^ = . Only the region below the red line is allowed. 


(4.15). Therefore, we need to take into account the window function Wr^ or the exponential 
suppression factor exp {—k'^/kjj) for large k modes. Since vt ~ kDiz 5 )~^, here, we simply 
take the Gaussian window function with r^, Wrj,{x) = {1 /{y/27rrT)^) exp(—r^/(2r|,)) and its 
Fourier component Wrj,{k) = exp(—r|,/c^/2), into account as a suppression factor for large k 
modes, which enables us to integrate these expressions explicitly, 


^AA,dc/g =1 \ ^ 0) 

25 


l-^/(/ + l) + 0^| 


~10 




(27r2 


where we have used 


j dkk‘^jf{kqo)exp(^-K^k‘^'^ =| 


TT 


2 ( .2 ) (^^2 

v '/o 


2V2qlrT V 

and the asymptotic expression for the modified Bessel function given by 

l{l + 1 ) 




V2 


-KZ 


1 - 


2z 


+ 0{z 




(5.7) 


(5.8) 


(5.9) 


Note that the spectral indexes are 4, which implies that the dimensionful power spectra 
P^^^{k —7- 0) for small k modes have no k dependence. We have also substituted rx — 
C)(0.1Mpc). 
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Figure 2. The figure shows the constraints on B-ti = /nl'Pk(^o) and trispectrum T-ji = ■ 

The contours of (2C'^'^ + C'^'^)/C'^® are inserted at the interval of 10 in the case with log;^Q ki/kf = 4, 
and a red dashed line is the line along which the non-linear correction is equal to Only the 

region below the red line is allowed. 


In passing, one also finds that (4.10) and (5.7) yield the weak constraints on power 
spectrum at /c ~ koiz^) 

iP7e(O(10Mpc-i)) < 0(0.1). (5.10) 

However, we have already known the stronger constraints on the scale by spectral distor¬ 
tions [7-13]. On the other hand, the constraints on g'NL can be about 10^ times weaker than 
those for tnlj comparing the cases with scale-independent local type tnl and c/nl- 

6 Conclusions and discussions 

We have discussed the inhomogeneities of the acoustic reheating developing the previous 
works which have mainly focused on the homogeneous part. Such inhomogeneities arise 
from the higher order correlation functions of primordial curvature perturbations because 
the acoustic reheating is a non-linear phenomenon, which first appear at the second order 
in perturbation. Produced (secondary) temperatnre perturbations are, in general, indis¬ 
tinguishable from the standard temperature perturbations linearly coming from primordial 
curvature perturbations. Given the fact that the linear theory prediction well matches with 
the CMB observations, the non-linear corrections should be subdominant compared to the 
observational i.e. the order of at large scales. Based on this condition, we 

gave constraints on higher order correlation functions, especially, the scale dependent local 
type non-Gaussianities of primordial curvature perturbations, in which large- and small-scale 
perturbations are correlated. These configurations of bispectra of primordial curvatnre per¬ 
turbations are complementary to those directly inferred from the three point correlation 
functions of the GMB anisotropies. Though the power spectra of the CMB anisotropies are 
used to constrain the (small-scale) higher order correlation functions of primordial curvature 
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perturbations, the extension to the higher order correlation functions such as bispectrum and 
trispectrum of the CMB anisotropies to further constrain primordial curvature perturbations 
is straightforward. 

In this paper, we have used only the total temperature perturbations to constrain the 
small-scale non-Gaussianities. However, the acoustic reheating can generate isocurvarture 
perturbations because the ratio of the number density of photon and baryon (dark mat¬ 
ter, neutrino) is perturbed universally, as already pointed out in [19]. The constraint on 
isocurvature perturbations may give more stringent constraints on the small-scale primordial 
perturbations. This topic will be discussed elsewhere [39]. 
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A Scale-dependent non-Gaussianity 

In Fourier space, the bispectrum of the comoving curvature perturbation is written as 

(77fe^77fc2 77fc3) = {2-K)^5^^\ki +k2 +k^)BTi{ki,k2,kz)., (A.l) 

and the local type one is parameterized as follows. 

k 2 , ks) = -^/nl(^i, h, 2perms.), (A.2) 

where Xa = ka/kQ. Now, let us consider scale-dependent local type non-Gaussianity defined 
as [36] 

fSt{k,,k2, fca) = 

with Xg = (xiX 2X3)3 and Xa = (x 2 -|- X 2 -|- xs)/3. is the /nl parameter at the pivot 

scale ko- In our cases, we can pick up the configuration of squeezed isosceles triangles, k <^ki, 
so we obtain the following forms. 

n f 

f§Liki,ki,k) = , (A.4) 

/NL(fci,^i>^) -• (A.5) 

For the top hat type, using the top hat function dehned as W{x) = 9{x — Xi) — 9{x — Xf) 
with 9{x) being the step function, f^^^^^{ki,k 2 ,k 3 ) can be written as 

^tophat(^i^ A; 2 , ks) + f}ll{W{xi) + 2perms.) 

+ .^L {W{xi)W{x 2 ) + 2perms.) 

+ jflw{xi)Wix2)W{x3). (A.6) 
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Assuming that the squeezed isosceles triangle configurations are enhanced, that is W{x) ~ 0, 
we can simplify (A.6) to 

h, k) = + hLWixi), (A.7) 


where /nl = 2^^ + trispectrum is also introduced in the same manner with 

the bispectrum, and we do not repeat it here [37]. Assuming the power law type, the scale 
dependence can be written as 


rNL(a:i, a:2, x[,X2) = 


G/A’ 


(A.8) 


where we have defined Yq = {xiX 2 X'iX/C) 4 and Ya = (xi + X 2 + xs + X4)/4. For the geometric 
averaged type, the above formula can be reduced to 


^"^(^1) ^1) ^1) ^i) 


CMBnr/2 , nr 12 

^NL 


(A.9) 


Top hat type is also defined as, 

(1T(xi) + W{x'i)) +rNLlT(xi)fF(xi). (A.IO) 
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